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Abstract. We consider the Starobinsky inflation with a set of higher order corrections 
parametrised by two real coefficients Ai , A 2 . In the Einstein frame we have found a po¬ 
tential with the Starobinsky plateau, steep slope and possibly with an additional minimum, 
local maximum or a saddle point. We have identified three types of infiationary behaviour 
that may be generated in this model: i) inflation on the plateau, ii) at the local maximum 
(topological inflation), iii) at the saddle point. We have found limits on parameters A; and 
initial conditions at the Planck scale which enable successful inflation and disable eternal in¬ 
flation at the plateau. We have checked that the local minimum away from the GR vacuum 
is stable and that the field cannot leave it neither via quantum tunnelling nor via thermal 
corrections. 
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1 Introduction 

The Starobinsky inflation [1] is an f(R ) = R + R? /6M 1 2 * theory and together with the original 
Guth model is the first model of cosmic inflation. It can be expressed as a Brans-Dicke theory 
with the Jordan frame potential proportional to (ip — l) 2 . After 35 years from its creation it 
is consistent with the recent CMB data [2] and it has developed many generalisations [3-9]. 
Its predictions (up to the N~ 2 terms, where N is the number of e-folds until the end of in¬ 
flation) coincide with the so-called Higgs inflation, where inflation is generated by the scalar 
field with V oc \q i > 4 and with non-minimal coupling to gravity of the form £(j) 2 R. 

Recently it has been shown [10-14] that the Starobinsky model is a part of an attractor 
of results on the (n s ,r) plane of several inflationary models based on theories of modified 
gravity. Its potential (as well as the potential of the Higgs inflation) may be modified by the 
higher order corrections, which take the form of the additional higher powers of Starobinsky 
potential in the action [15, 16]. The higher order terms may generate a steep slope of the 
Einstein frame potential which limits the Starobinsky plateau and the amount of e-folds that 
can be generated during inflation. In this paper we want to investigate those generalisations, 
especially in the context of possible existence of local minima and maxima of the potential 
and their influence on the evolution of the inflation, existence of the successful inflation and 
late-time evolution of the Universe. 

The structure of this paper is as follows. In the Sec. 2 we analyse the evolution of inflaton 
and space time for the certain generalization of the Starobinsky potential. In the Sec. 3 we 
investigate several types of inflation that could be generated by considered potentials. In the 
Sec. 4 we include non-zero temperature corrections. In the Sec. 5 we discuss the possibility 
of quantum tunnelling from meta-stable vacua of the Einstein frame potential. Finally, we 
conclude in the Sec. 6. 

In the following we use the convention 87 tG = M~ 2 = 1, where M p ~ 2 x 10 18 GeU is the 
reduced Planck mass. 
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2 Generalisation of the Starobinsky inflation 


2.1 Jordan frame analysis 

Let us consider a Brans-Dicke theory in the flat FRW space-time with the metric tensor of 
the form ds 1 2 = — dt 2 + a(t) 2 {dx) 2 . Then the Jordan frame action is of the form 


S 


d A x^/\g\ 


<pR-^ ( p) 2 -u{<p) 


+ S n 


( 2 . 1 ) 


where uj = const and S m is the action of matter fields. From now on we will assume that 
u = 0, which makes the model equivalent to the f(R ) theory. Then, for the homogeneous 
field tp the field’s equation of motion and the first Friedmann equation become [3] 


ip + 3 Hip + —((pUp — 2U) — — (pm ~ 3Pm) 


3 H + 




2<p) 4 

Pm + 3 H(pm + Pm) = 0 
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<P 
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ip ip 


( 2 . 2 ) 

(2.3) 

(2.4) 


where U v := and pm and Pm are energy density and pressure of matter fields respectively 
1 . Let us note that for p = 1 one recovers the general relativity (GR). Thus the <p = 1 will 
be denoted as the GR vacuum. 


The Starobinsky inflation is a theory of cosmic inflation based on the f(R) = R + 
R 2 /6M 2 Lagrangian, which can be generalized to Brans-Dicke theory with general value of 
ui. The Jordan frame potential of the Starobinsky model is following 

U s = \m 2 {v-1)\ (2.5) 

where M is a mass parameter, whose value comes from the normalisation of the primordial 
inhomogeneities. For co = 0 one finds M ~ 1.5 x 10 -5 . The Starobinsky potential is presented 
at the Fig. 1. In this paper we consider the extension of this model motivated by the Ref. 
[15], namely 

u = Us i 1+Xl w* +X2 Mj>) • < 2 - 6 > 

where Ai, A 2 are numerical coefficients. In order to avoid U —>• —00 for ip —> 00 we assume 
that A 2 > 0. The sign of Ai is undetermined. The motivation for the existence of the higher 
order terms is as follows: we assume that the Jordan frame field is a singlet of some theory 
beyond the Planck scale. Integrating out the heavy degrees of freedom gives M p in the 
denominator of the higher order terms in the potential. In addition, we require that the GR 
vacuum does exist in order to restore classical gravity. 

The model as an f(R) theory 

As mentioned previously one can define the f(R) function (where R is the Ricci scalar) for 
which 

ip = F{R):= dR ’ U{<p) = \{RF-f), U V = R . (2.7) 

1 In this paper we refer as matter fields to all perfect fluid components of the energy-stress tensor, like dust, 

radiation or scalar fields. 
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Thus, from p(R) one can reconstruct f(R) by / = f p(R)dR. For U = Us one finds 
f{R) = R + For the general values of A 2 one cannot solve U v = R/2 to obtain 

( p = p(R) = F(R), because there is no general solution of the fifth order equation. In such a 
case we cant write an explicit form of f(R). On the other hand, for A 2 = 0 and Ai = A one 
finds 

**> = *■(*)-I + < 2 - 8 > 

f(R) = (4A 2 (1 + 12RA) - (B - 6RX 2 ) A - (M 2 X + 3 RB) A 2 ) , (2.9) 

where A = (6RX 2 + 2^X 3 (2Af 2 + 9i? 2 A)) M“ 4 / 3 and B = sjX 3 (2Af 2 + 9R 2 X) - 3RX 2 . 

For A —> 0 one recovers the Starobinsky model. Considering a general value of u ; shall be the 
next step in our future analysis. 

2.2 Einstein frame analysis 

The gravitational part of the action may obtain its canonical (minimally coupled to p) form 
after transformation to the Einstein frame. Let us assume that p > 0. Then for the Einstein 
frame metric tensor 

9nv = V9nv , dt = s/pdt , a = ^Jpa (2.10) 

one obtains the action of the form of 

T Sm[gfiui A\ 1 (2-11) 


S[g^,<p\ = j d A x^-H 


Io_3 
2 4 1 p 


EM 


where V is the derivative with respect to the Einstein frame coordinates. Matter fields are 
now explicitly coupled to p due to the fact that d A Xy/^g = p~ 2 d A x^U-g. In order to obtain 
the canonical kinetic term for p let us use the Einstein frame scalar field (j) 


<t> = 




( 2 . 12 ) 


The GR vacuum appears at cj> = 0. The action in terms of g^ and 0 looks as follows 


S= d A x-\f--fj 


-R — - ( V 
2 2 


<py ~v{ 4 >) 


T 5 'm \_9fiv > tp > 


(2.13) 


where 


V = 


U{p) 


p+ 


ip=<p(<P) 


(2.14) 


and R is the Ricci scalar of g In this section we will assume that the space-time may be 
described by the flat FRW metric tensor. Let us define the Einstein frame Hubble parameter 
as 


ft , 

a 


where 


a : = 


dd 

dt 


(2.15) 
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Then for pm = Pm = 0 the first Friedmann equation and the equation of motion of (j) are 
following 

3 P 2 = ^ ,2 + V^), (2.16) 

4>" + m<\J + v 4> = o 1 (2.i7) 

where V* = %■ 

Minima and maxima of the Einstein frame potential 

The Einstein frame potential obtained from the Eq. (2.6) has following features: For 
Ai, A 2 > 0 one obtains the Starobinsky plateau, which ends with a steep slope when one 
of the higher order terms starts to dominate. In such a case V has no maximum and the 
only minimum appears at 0 = 0. Thus cj) ends its evolution at cj) = 0 for any set of initial 
conditions. For Ai < 0 the potential has a minimum at <f> = 0 and may have a local maximum 
at 4> max > 0 and a minimum between the plateau and the slope. The last minimum will be 
denoted as d> min (where 4>mm > 4> max) and its depth will depend on the relation between Ai 
and A 2 , so </> m i n can be a local or global minimum. Potentials with minimum at (j) = (f > m ; n are 
presented at Fig. 1 and 2 

Let us assume that Ai < 0 and A 2 > 0. In order to find an extreme value of V one needs 
to solve the equation = 0. For <f> of order of a few M p (which is the case for 4> m \ n or r/> max ) 
one finds exp(y / 2/3</>) 3> 1, so the condition = 0 simplifies into 

8 + 6e 3 V^M 2 Ai + 9e 5 V^M 4 A 2 = 0 . (2.18) 

This equation may have solutions as long as the Ai can dominate other terms for some values 
of (j). In such a case, for Ai < 0, the minimum (maximum) or the saddle point exists for 
< 0 (Vfptp > 0) or V^j, = 0 respectively. Let us consider the latter case in which has 
a saddle point at some (j) = (f> s , but does not have any maxima or minima besides the one 
in (j> = 0. Then, assuming exp( y/2/3 4>) 1 and Ai M 2 A 2 one finds that the condition 

= 0 simplifies into 


- 2 + 3AiM 2 e 3 


+ 9A 2 M A e 


= 0 


(2.19) 


at (j) = (j) s . Then from Eq. (2.18,2.19) one finds 


4*s 



3Ai M 2 ) ’ 


5A 3/5 M 2 / 5 
- 2V532/5 ’ 


V. = \M* b - P (24M%) ‘/A , 


( 2 . 20 ) 

where \ s fixes the relation between Ai and A 2 in order to provide the existence of a saddle 
point and V s = V((j) s ). Let us note that for <f> ~ (f> s one obtains a saddle point inflation, 
which in principle could significantly decrease the scale of inflation. The for Ai = — A s is 
presented at the Fig. 2. For |Ai| < X s the Ai term is always subdominant compared to the 
other terms in V(4>) and it can be neglected in the analysis. For Ai > — X s the potential has 
no stationary points besides the minimum at 4> = 0. Lets us assume that Ai < — X s , which 


- 4 - 




means that 4> m i n and L max do exist. At the minimum exp(y / 2/30 m ; n ) ~ 0(M x ), so the 
first term in the Eq. (2.18) is negligible. Thus, for |Ai|, A 2 < M^ 1 one finds 


- o\/o !og -o 


2 Ai 


3 M 2 X-2 


Assuming Ai, A 2 < M 1 one finds the following value of V m 


Enin : — V (</> m in) — 


3 M 2 (4A 2 - A?) 
I6A2 


( 2 . 21 ) 


( 2 . 22 ) 


Thus E m i n > 0 Ai G [—2^/A^, 0). If the field oscillates around d> min it has a mass defined 
by 

A^ Af2 

ml := 1440 = d>min) -~7— • (2.23) 

M 

One can apply the same procedure to calculate </> max . At the plateau 1 <C exp(y / 2/3 fa) <C 
M^ 1 , which makes the third term of the Eq. (2.18) negligible. Thus <^> max and E max := 
E (</>m a x) 



-4 


V 3M 2 Ai 



(2.24) 


Let us stress once again that 0 max is only a local maximum of the potential. The values of 
d> min and cj > ma x are plotted at the Fig. 4. One could argue that typical values of Ai and A 2 
are of order of unity, but in general they may be much bigger. The only limitation is that 
the steep slope shall begin for </> > 5.5, so the plateau is long enough to generate successful 
inflation. As we will show this requirement gives |Ai| < M _1 and A 2 < M~ 2 . 


Evolution of 4 > 

The evolution of the field is as follows: the field starts its evolution at fa > ^> m in and rolls 
down towards 0 m j n . If its kinetic energy is sufficiently big it rolls up to the plateau and gen¬ 
erates Starobinsky-like inflation. The length of the plateau is typically around fa — <pf ~ 10, 
where fa and <f>f are the beginning and the end of the plateau respectively. This is more than 
enough, since in Starobinsky inflation one needs 4>i — 4.7 in order to generate 60 e-folds of 
inflation. Another interesting issue is that around the middle of the plateau the potential 
starts to decrease with 0, due to the existence of local maximum at </> ma x- As we will show 
this could lead to eternal inflation [17] at d w x- To avoid that one could require that the 
plateau between <f> = 0 and <^ max is long enough to generate successful inflation. This require¬ 
ment is fulfilled for all |Ai|,A 2 ~ 0(1). Nevertheless for |Ai|, A 2 3> 1 one can obtain a short 
plateau for which the last 60 e-folds of inflation are partially generated by the inflaton in the 
0 > regime. At the end the field rolls down to cf) = 0, oscillates around the minimum 
and reheats the universe. 


Let us assume that the field is somehow bigger than </> m i n and the A 2 term dominates 
the potential. Then E ~ A 2 exp(—2^2/3 fa)Ug oc exp(4y / 2/3 fa). Such a potential cannot 
generate the slow-roll evolution of the held since the slow-roll parameter e = 16/3i$>l. If 
the held cannot escape the minimum at </> m i n it oscillates around it and due to the existence 
of the cosmic friction hnishes its evolution at <f> = 0 m i n . If E(0 m i n ) > 0 then one obtains 
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an exact de-Sitter solution for the evolution of the space-time. This kind of inflation suffers 
from lack of consistency with PLANCK [18] and BICEP [19] data because of the perfectly 
flat power spectrum of primordial inhomogeneities. Another problem is a lack of the graceful 
exit, which ends inflation and allows to reheat the universe after inflation. A possible solution 
to that problem may be the quantum tunnelling, which will be discussed in the following 
parts of this paper. The possible solution to the flat power spectrum problem could be to 
introduce the additional curvaton field responsible for the generation of the primordial inho¬ 
mogeneities. The potential of the curvaton should not be perfectly flat in order to generate 
perturbations consistent with observational data. Thus, the de Sitter inflation should not 
last for too long. Otherwise the curvaton would vanish before the tunnelling. Another issue 
appears when (fi m i n is a global minimum of V. Then, the field which managed to reach (f >o 
may tunnel to 4> m i n , which changes the effective value of gravitational constant. In such a 
case the theory would not recover GR at low energies. 

Fig. 3 shows the areas in the phase space of initial conditions that lead to <ft —> <t > min at 
the late times. The region with enables <f> to reach </> = 0 lies in the 4>o K(^>o) regime, so 
the kinetic term domination is required to reach the GR minimum via the classical evolution 
of the field. Let us note that no positive value of <f >q allows <ft to reach the GR minimum. 




Figure 1. The Einstein frame Starobinsky potential (thick black line) and its generalisation from 
the Eq. (2.6) for M = 1.5 x 10~ 5 and for different values of Ai and A 2 . Left panel: The potential 
with minimum, saddle point and without stationary points for (f> > 0. Right panel: The depth of the 
potential for different values of Ai. Local minima with Enin > 0 may generate the de Sitter solution 
without graceful exit. 


Planck scale limit 

One can ask when the initial conditions give the potential energy of order of the Planck scale. 
For given values of Ai and A 2 let us define cf> p i such as V(4> p i) = M Whenever 1 f> > cj) p i the 
potential term is bigger than the Planck scale and quantum gravity corrections should be 
included. The potential V(4>) reaches the Planck limit when the A 2 term dominates. Thus, 
one can easily assume that in that limit U ~ A 2 Ug. Then, from the definition of cj) p i one finds 

= \/l log ■ (225) 
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Figure 2. Left Panel: For fixed relation between Ai and V ^2 one finds the same value of F m i n for 
all A 2 . Right panel: for Ai = — A s . We have assumed A 2 = 10 fc , where k £ {— 7, — 6,..., 3,4} 

(smallest/biggest A 2 are the right/left orange curves). Black dots represent <p = <f> s . Both X s and (j) s 
are taken from Eg. (2.20) and they appear to be very accurate. One finds V^cfs) ~ A^M 18 / 5 1. 


The biggest contribution to <f> p i comes from M. On the other hand the A 2 dependence is very 
weak. For all realistic values of A 2 and for M = 1.5 x 10 -5 one finds 4> p i — 20. 

The Planck scale as the limit of the classical theory is not set precisely. Different models 
of quantum gravity predict different scales of validity of the classical theory. For instance in 
loop quantum cosmology (LQC) [20], the maximal energy density of the system, at which the 
quantum corrections are the most significant is of order of 250M p . This would allow bigger 
values of 0, which lie within the classical limit of the theory. The situation becomes more 
complicated when one considers LQC corrections to modified theories of gravity [21, 22], 
Then, the maximal energy scale depends on the frame in which the system is quantized, and 
on the value of 4> at the moment of maximized energy density. Thus, one cannot set one 
scale that would clearly limit the 4> from the point of view of the validity of the classical 
theory. The <f> p i is just a realistic approximation, especially since the slope is very steep and 
any increase of <f> over <f> p i increases the energy density of the system by orders of magnitude. 
The issue of the LQC corrections to this model shall be a future extension of the analysis 
performed in this paper. 

The problem of initial conditions 

Let us note that due to stronger limits on the tensor-to-scalar ratio r the PLANCK data 
favour the plateau-like potentials (like in Starobinsky inflation) and disfavour the power-law 
potentials (like m 2 (j) 2 ). Nevertheless, the plateau-like models of inflation struggle with some 
difficulties mentioned in the Ref. [23] . One of the problems, which appears within the anal¬ 
ysis of the Starobinsky model, is the fact that the Einstein frame potential is limited from 
above by the scale of the order of M 2 <C M p . Thus, if one would set initial conditions at 
the Planck scale the potential term would be always subdominant and the (9*0) 2 term may 
dominate the universe, which would lead to strong inhomogeneities. While potentials with¬ 
out the upper bound generate the Planck scale initial conditions (<9j0) 2 ~ <f> 2 ~ V(f >) 2 , for 
which the potential term can dominate the evolution of the field and homogenize the universe. 


2 For argument see e.g. the Ref. [24]. 
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Figure 3. Both panels present regions in the phase space of initial conditions, for which the GR 
minimum cannot be reached. Red curves correspond to po = M£ and po = 100 (dashed and 
dotted red lines respectively). As shown in the right panel the field can reach the GR vacuum only for 
</>q 2 V((j)tf) and only for <f >' 0 < 0. Even super-Planckian values of <f' 0 >0 would give the evolution 
which ends at <p m i n . 



Ai 


Figure 4. The plot show numerical results for (f> m in and (/> max (solid and dashed lines respectively) as 
a function of Ai for several different values 0 /A 2 . Not surprisingly <f max does not depend on A 2 , which 
dominates the potential for bigger values of <f>. Both 4> m , n and 4> max meet at <f s when the |Ai| > A s 
is violated. Those results are perfectly consistent with analytical analysis of </> m j n and </> max (see Eq. 
(2.21,2.24)) besides the narrow region of very small Ai, where the A 2 term shall also be taken into 
account in the Eq. (2.24). 


For the Starobinsky potential with higher order corrections this problem seems to be 
weakened. The steep, exponential slope at big <f> reaches the Planck scale around <j) ~ 20 and 
the potential can have significant contribution to the initial energy density. Nevertheless the 
problem of initial conditions is still an issue of our model, since inflation cannot start on the 
slope and inhomogeneities can still dominate the evolution of the universe when the field is 
on its way to the Starobinsky plateau. The other issue appears for Ai < — A s . In order to 
pass local minimum and reach the plateau we require the kinetic term domination at the 
Planck scale. Thus, before reaching the beginning of inflation the inhomogeneous terms may 
become significant and spoil our analysis based on FRW metric. 
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Let us note that the slight modification of the Starobinsky potential based on the 
f(R) = R + aR n , where n < 2 is even better cure for the problem of initial conditions. For 
n/2 the Einstein frame potential is of the form [ 6 , 7] 


V oc Vs x 





4>»l 


exp 


2 — n 
n — 1 



(2.26) 


where Vs := Us/'p 2 - This potential is not limited from above, so V(cj)pi) can have large 
contribution to the energy density at the Planck scale. Another advantage of this model 
is that inflation can start for any cj) and therefore to obtain successful inflations without 
significant inhomogeneities one needs only one bubble at the Planck scale in which V(<p) 
dominates. 


3 Inflationary dynamics 


Inflation with and beyond the slow-roll approximation 

In this section we will discuss the slow-roll approximation and inflation in the Einstein frame. 
Inflation takes place on the plateau, where the influence of the A 2 term is very small, but as 
we will show there are also other mechanisms generating inflation, which strongly rely on A 2 - 
We want to investigate how non-zero values of A parameters deviate the inflation away from 
the Starobinsky model. Let us assume that 4>" -C V^. Then, for pm = Pm = 0 one obtains 

+ v+ ~ 0, 3 n 2 ~v. (3.1) 


This approximation holds for non-zero values of V^, so one cannot use Eq. (3.1) at <f> = 0 max . 
Thus we will use the slow-roll equations for cj) 7 ^ ^max- The cosmic inflation takes place as 
long as following slow-roll parameters are much smaller than one 


e := 


1 

2 




(3.2) 


The number of e-folds generated during the inflation is in the slow-roll approximation equal 



where indexes i and / refer to initial and final moments of inflation respectively. Namely, 
ti is the first moment when both slow-roll parameters are smaller than one and tf is the 
moment when any of slow-roll parameters becomes bigger than one. 0 / is independent of 
higher order correction to the potential, but in general may depend on Ai, A 2 and initial 
conditions chosen for the evolution. 


For Ai < —A s , A 2 > 0 there is a local maximum at </> max , which may influence markedly 
the total number of e-folds generated during inflation. One can fine-tune initial conditions 
to fit the following evolution of the inflaton: </> has just enough energy to climb to the 0 max , 
so at the local maximum its kinetic energy is almost zero. Then cj) stays very long around 
4> m ax , which in principle may generate unlimited number of e-folds according to the classical 
evolution of 0. Even if Ai, A 2 are so big that the Starobinsky plateau is too short to generate 
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at least 60 e-folds of inflation, one could generate sufficiently large N during the evolution 
around cp max . The infinite production of e-folds appears for fine-tuned initial conditions which 
lead to <p = cp ma x and cp' = 0. 

This effect is limited by the existence of quantum fluctuations of the inflaton. The 
vev of the inflaton evolves due to the classical evolution and due to quantum fluctuations. 
During one Hubble time those effects modify the value of the field by 6 (pc ~ ( P'/H and 
Wq 'H/2'k respectively. Whenever 5cpQ > 5cpc the classical equations of motion cannot 
be used to describe the evolution of cp. In our case this means that whenever the held 
rolls up or down on the attractor trajectory around cp max and satisfies ScpQ > 5<pc its vev 
is pushed after one Hubble time towards (pram or GR vacuum by quantum fluctuations of 
the inflaton. To see how efficient quantum fluctuations are let us define the region on the 
attractor which gives 5c/)q > 5<pc by cp G (<p max — Acp, cp m a.x+Acp), where Acp satisfies condition 
5cpQ (^max =t A cp) = 5<pc{(P m a .x ± Acp). Since the Eq. (3.1) is satisfied at the attractor one 
finds _ 

=* A(/>~ ^-(-M/6Ai) 1/3 . (3.4) 

For Ai = —1, M = 1.5 x 10 -5 this gives Acp ~ 5.4 x 10“ 4 , which is two orders of magnitude 
bigger than the typical value of 5<Pq for this particular value of M. During one Hubble time 
from our comoving Hubble radius one obtains e 3 causally independent regions. In some of 
them dcpQ < 0 and the held is pushed toward the Starobinsky plateau. Nevertheless one 
would need 100 such a fluctuations (or few with extraordinary amplitude) to leave the part 
of the attractor with quantum huctuations domination. It is indeed highly improbable event 
and in most of horizons generated while | <p — (/> max | < Acp it shall never happen. Thus, 
whenever the held reaches the attractor at | cp — <^> max | < Acp the held could never reach the 
classical evolution domination regime. Thus one obtains the eternal inflation. To avoid this 
problem one would need Ai > 3 x 10 -4 M -2 ~ 10 6 , which gives Acp < 5<Pq. In such a case 
cp would leave the regime of quantum huctuations domination after t ~ 1/H. The other 
issue is that for |Ai|, A 2 —> 0 one hnds that Acp and cp m ax are growing like A 1 7 and — logAi 
respectively. Since power-law functions are growing faster that the logarithm one expects the 
whole plateau to be inside the (cp m ax — Acp, cp ma . x + Acp) region. Fortunately this behaviour 
of Acp does not hold for |Ai|, A 2 <C 1 and Acp region is just swallowed by the other region of 
eternal inflation at cp > 17.3. In principle the quantum huctuations around the maximum 
could lead to the so-called topological inflation. This issue was partially analysed in the Ref. 
[16] and will be extended in our further work. 

As shown at the Fig. 5 the inflaton during its evolution reaches the attractor separated 
into two regions. The one that exists for cp < cp max leads to Starobinsky-like inflation and 
to oscillations around the GR vacuum. Initial conditions, for which cp reaches the attractor 
at cp > </>max give the evolution which always ends at <p = cp m j n and cp’ = 0. Those regions 
of the attractor are separated by the region where quantum huctuations dominate around 
0maxj where classical equations of motion do not hold. 

As mentioned previously in this section whenever the inflaton reaches the attractor 
at \cp — cp mstx \ < Acp one obtains the quantum huctuations domination regime and eternal 
inhation. To avoid that let us assume that cpi < cp mSLX — Acp. If cpi < log (—M _2 Aj" 1 ) /\/6 
then one can expand the number of e-folds in terms of the mass M. Then, for cpi = cp one 
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Figure 5. Both panels present the phase space analysis for Ai = — = —1. Different initial 

conditions give the evolution, which leads to one of two regions of the attractor: a) The blue one gives 
inflationary trajectory at the Starobinsky plateau and oscillations around (f> = 0. b) The red one leads 
to </> m in and the De Sitter solution. Dashed purple line represents the slow-roll approximation, which 
describes the evolution the attractor from <f> ~ 1 to </> m in- Axes originates at (0 m ax>O). 


obtains 

N ~ ^ ^3eV^ - y/6(/> - ^e 4 V^M 2 A 1^ . (3.5) 

In fact the biggest amount of e-folds can be generated when < f> is closer to c/> max — A0 since 
in this regime is smallest. In such a case the linear approximation collapses and one need 
to employ numerical methods. Keeping that in mind let us define the maximal amount of 
e-folds generated during the classical evolution of the field by 

/'^max T/~ 

N max = / tt#. (3.6) 

J<h v <t> 

Its values for different Ai are presented at the Fig. 6. We require that IV max > 60, which 
gives Ai < 5 x 10 5 . In fact Ai could be much bigger but this would require Ai = —A<j and a 
saddle point inflation. One can see that for Ai < —A s (right from the peak for given A 2 at 
the Fig. 6) the lV max does not depend on A 2 - This comes from the fact that 0 max and A <f>, 
which determinate maximal allowed 0j, depend mostly on Ai and very weekly on A 2 - 

The primordial inhomogeneities generated in this model do not deviate much from the 
Starobinsky inflation as long as </> max is separated from the scale of freeze-out of perturbations 
by few M p . In such a case the last 60 e-folds of inflation happen in the regime of negligible 
influence of higher order terms. On the other hand, for <f> mSLX ~ 5 or for the saddle point 
inflation one expects significant changes in the shape of power spectrum and in the value of 
M as compared to the Starobinsky model. This issue is currently under investigation and it 
shall be discussed in a forthcoming paper. 

Eternal inflation for large <j> 

As mentioned above the evolution of cj> may, under some circumstances, be strongly influenced 
by quantum fluctuations. For the Starobinsky inflation one finds that 

5c/)q > 54> c for 4> > (f c = -/flog 
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Figure 6. Left panel: Number of e-folds generated during inflation as a function of initial 
field value. The Ai = A 2 = 0 corresponds to Strobinsky , and other lines illustrate how fast 
is the departure from this simplest case due to corrections. The value of A 2 in all cases with 
negative X± was fixed be requiring that the value of the potential in the minimum is M 2 / 10. 
Right panel: iV max for different values of Ai and A 2 , where Ai < 0. Peaks correspond to the sad¬ 
dle point of the potential at Ai = — A s . Flat lines left from peaks: minimum and maximum do not 
exist and inflation may start as soon as the slope ends on the plateau. The A 2 = 1 line right from 
peaks: inflation may start as soon as <f> < </> max — A0. All the lines converge to the same function for 
large enough Ai because after the new minimum appears number of e-folds to the leading order is only 
a function of Ai 


For M = 1.5 x 10~ 5 one finds 4> c ~ 17.3. For the Starobinsky inflation with > <f> c one 
obtains the uplifting of <j) by quantum fluctuations, which disables the field to roll down. In 
such a case inflation never ends, which is the so-called eternal inflation. As mentioned in [15] 
the higher order terms in the potential may solve this problem. Indeed, if any of higher order 
terms would generate the slope or the minimum at some (f> < (f> c one would be free from the 
eternal inflation. To satisfy this condition one needs 


Ai > —^—e 2 \/~3 ^ c — 4 ~ 3 2 x 10 -3 

Al> 3A P 6 ~32^- 3 - 2XW 

This result does not depend on M. 
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B -V§*c = 2~ w tt- 4 


A2 > 9W 6 ’ V ^ = 2 ~ iU7r ~" ~ 2 ‘ 4 X 10 


-4 

(3.7) 


4 Thermal corrections to the evolution of the system 
Introduction 

When the field is trapped in d> min with V(<f> min ) = 0 one can try to evolve it to (f> = 0 with the 
help of the thermal corrections. The non-zero temperature of the universe may be obtained 
by the dissipation of the inflaton’s energy density into relativistic degrees of freedom from the 
standard model (SM). The idea of thermal corrections to the potential has been developed 
in the context of the “Mexican hat” potential [25, 26], where structure of the vacuum breaks 
spontaneously the symmetry of potential. While the temperature it growing the potential is 
modified due to the existence of the effective potential term proportional to T 2 ^ 2 . This ad¬ 
ditional term fills the minimum, pushes the field towards zero veg and restores the symmetry 
of the potential. 
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In our case there may be two sources of non-zero temperature in the inflationary Uni¬ 
verse before the reheating. First of all one can reasonably assume that the universe before 
inflation was warm or that the inflaton has decayed into relativistic degrees of freedom. The 
decay of the inflaton at </> = (/> m i n is in principle the same as the post-inflationary reheating 
of the universe. The <j) is coupled to scalars, fermions and vectors of the SM (or to other 
particles, which can decay into the SM degrees of freedom), which are being produced during 
the oscillations of the field around d> min - If V(d> min ) > 0 then the vacuum energy generates 
the de Sitter space-time, in which the energy density of relativistic particles will decrease like 
e -4Ht . Thus, the universe would remain cold and thermal corrections will not modify the 
potential. 

Let us define the Einstein frame potential with non-zero temperature by V(</>, T) = 
V(cf>) + AU(i tf>,T), where A V(cj),T) is the thermal correction to the potential. The thermal 
correction to the Einstein frame potential has a following form [27] 

Av= S v «, t 4 - 1 ) 

where T is the Einstein frame temperature. This means, that the Einstein frame radiation 
energy density is equal to pn = 7t 2 /30 gvT 4 , where g * is the number of relativistic degrees of 
freedom. In general 5 * = g*(T), since g* may decrease with the energy scale. 

The features of the Einstein frame potential with the thermal correction 

In this section we will apply the procedure introduced in the Sec. 2.2 into the V(</>, T) case. 
As previously, we will use the fact that exp(y / 2/3 ~ M ^ 1 and we will assume that 

| Ai|, A 2 < M~ l . Since pr <C M 4 we also assume that T < 1. Then V has a minimum at 
</>min defined by 



2Ai \ 

3A 2 M 2 ) 



brain ~ ~ \ I ~T Z 


(4.2) 


One can see that the correction is rather small and it does not significantly shifts the mini¬ 
mum. The potential energy density at <j) = </>^ in is following 


V„ 


3 M 2 . M 2 A? o 

(4A 2 — A x ) + ^ ,— T — Vmm 


16A 2 


24A 2 


+ 


M 2 A 2 r2 
24A 2 


(4.3) 


The non-zero temperature is uplifting the potential in the minim. However, this effect is way 
to small to fill the minimum and to push (j) towards the GR vacuum. Even if one would lift T 
towards the Planck scale the minimum at wouldn’t be filled. Higher order corrections in 
T would enter the Eq. (4.3) and the energy density would start to decrease with temperature. 
This issue is depicted in Fig. 7. 


The thermal corrections could be also treated as a way to make the minimum at 4> m i n 
more shallow while the field is on the slope. In such a case the region in the phase space of 
initial conditions which allowed (j> to reach the GR vacuum would be bigger. Unfortunately, 
due to the fact that AU <C V for T <C M p , this effect would be very small. Thus, the non-zero 
temperature of the universe shall not significantly change the evolution of the system. 
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Figure 7. Both panels present V(<j>,T) for Ai = —1, A 2 = 1 and for different values ofT. Left Panel: 
maximal realistic values of T. The minimum is being uplifted, but the inflaton still cannot escape 
to the GR vacuum. Right panel: even for unrealistic values of T the minimum will not be filled by 
the thermal correction. For T > 2>M p the V m i n starts to decrease towards the negative values. These 
results remain the same even for |Ai|, A 2 1. 


5 Quantum tunnelling from meta-stable minima 


Whenever a new minimum appears in the potential due to higher order corrections (i.e. for 
Ai < —A s ), the field may be unable to cross the barrier due to classical evolution if its initial 
speed is to small. In this section we will discuss whether it is possible for the field to cross 
the barrier due to quantum fluctuations. 

We will focus mainly on two limiting cases for such an event. First one is the Hawking- 
Moss effect (HM) [28] describing a simple fluctuation of the field, large enough for it to pass 
the barrier. The second case is quantum tunnelling. Here we will use the standard formal¬ 
ism of Coleman and De Luccia (CDL) [29], assuming that vacuum decay proceeds through 
nucleation of true vacuum bubbles within our false vacuum. Details of this calculation are 
described in the following subsection5.1. During this analysis we show two examples of po¬ 
tentials with fixed values of the false vacuum minimum H m ; n = M 2 /10 and V m \n = M 2 /100. 
In the remaining part of this section we use these conditions and calculate A 2 for each value 
of Ai. While we use a precise numerical result, the analytical estimate coming from (2.22) 
will be sufficient to understand these results, namely 


A2 


A? 


4 M _y . 

^ 3 W 1 


(5.1) 


Existence of the plateau between the two vacua is crucial in vacuum stability analysis 
because CDL instantons cannot penetrate a barrier which is too flat. The crucial parameter 
here is [30] 

O Q ^(^(^ max) 

^ — V 

v max 

for f3 > 4 both CDL and HM solutions exist. We can immediately see a problem because this 
parameter is proportional to rj slow-roll parameter from (3.2). This suggests we can either 
have CDL solutions or large number of e-folds during inflation, not both. This suspicion is 
confirmed by Figure 8 which shows parameter /3 for our two example scenarios. 


(5.2) 
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Figure 8. Parameter f3 from (5.2) in our two examples. The CDL solutions exist only if j3 > 4. And 
below /3 = 1, only HM solutions exist. The value of |Ai| = 10 6 is the largest allowing 60 e-folds of 
inflation on the plateau. For such a big Xi the topological inflation at <fi = </> max is negligible, since 
Ac j> < H/2 -k. 


The barrier penetration through CDL instantons becomes possible only for very large 
values of parameters A. This happens because the larger the couplings are the shorter and 
less flat the potential barrier becomes. One can see that when ft reaches 4 these potentials 
are already excluded by lack of suitable inflation as discussed in Section 3. The only way 
to generate over 60 e-folds for such a big A parameters would be to assume that Ai = — X s . 
However, such a potential would not have a minimum and thus the issue of quantum tun¬ 
nelling would not exist. 


For both type of instantons considered in this paper one needs to calculate the decay 
probability of the metastable vacuum defined by 

T = Ae~ s , (5.3) 


where S is an action of the considered type of instantons. The prefactor A is derived from 
quantum corrections to the bounce and always is less important than the exponent. Action 
of an HM instanton is simply the difference between the action of homogeneous solution for 
field in the false vacuum and on top of the barrier at </> max 


<ShM = 'S'l^’max] - S^min] 


-24?r 2 ( 24vr 2 
~V ^ V ■ 

v max v mm 


(5.4) 


Action of an CDL instanton is the difference between the instanton solution and the homo¬ 
geneous background (false vacuum), 


<ShM = S^CDL] - 5'[</ > m in] = S^cdl] + 


24vr 2 

Fmin 


(5.5) 


We will describe calculation of S^cdl] hi detail in the following subsection 5.1. However it 
is already clear that a CDL action is always smaller than the HM one because the instanton 
solution </>cdl interpolates between the two vacua rather than constantly taking the highest 
value on top of the barrier as in HM solution. 
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Figure 9. Action of CDL and HM instantons with the false vacuum fixed at V m i n = Al 2 / 10 (left 
panel) and V min = M 2 /100 (right panel). 

Figure 9 shows actions of both solutions in our example potential. The resulting lifetime 
of the vacuum can be bounded from below by choosing the Planck scale as a dimensionfull 
quantity describing our potential, 

r = -p—rjje. S ~ 1.8 x 10 25 e s s, (5.6) 

4v tM£ v ’ 

the current age of the universe is Tjj ss 4.3 x 10 17 s. Our result is huge compared to the 
current age of the universe despite Planck scale suppression, leaving absolutely no hope for 
a correct power spectrum after the tunneling finally occurs (as described in Section2.2). 

The only remaining loophole are oscillating bounces described in [30]. These bounces 
exist for values of (3 between 1 and 4, and so they may dominate vacuum decay when CDL 
instantons do not exist. However when both CDL and HM solutions exist, actions of oscil¬ 
lating bounces are between the two. Thus looking at Figure 9, one inunediatly sees from 
continuity of the solutions that the oscillating bounces are also completely irrelevant in our 
problem. 

In summary we conclude that in the described potential the metastable vacuum is always 
significantly to long lived. And all solutions involving field being trapped in such a minimum 
due to classical evolution are excluded. 

5.1 CDL instantons calculation 

We are interested in an 0(4) symmetric scalar field configuration </> = </>(t), with the Einstein 
frame metric given by ds 2 = dr 2 + r(r) 2 (dD) 2 . Here dCl is an infinitesimal element of the 3D 
sphere, and r(r) is the radius of that sphere. The resulting metric tensor is of the form of 
the FRW metric with the curvature parameter k = +1. Euclidean action in Einstein frame 
takes the form 

S E = 2vr 2 J drr 3 Q^> 2 + V + , (5.7) 

where R = 6 + (yr)^ — an d 4> r = ^- The equation of motion of the scalar field 

reads 

(j) TT + 3— (j) T = , (5.8) 

r 
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while the first Friedmann equation reads 


r T =^l + j ( 5 - 9 ) 

One can show that scale factor r crosses zero twice [31]. Without loss of generality we chose 
value of r of the first one to be r = 0, the other one we call r erK i. The appropriate boundary 
conditions then are 


•MO) = Mrend) = 0 (5.10) 

r{ o) = r(Tend) = 0 

The final initial condition needed for our equations is the initial field value fo. For initial 
values (j )o to close the minimum of the potential the held will pass the other minimum, while 
for values too far from it it will start oscillating near maximum of the potential by the time 
when r crosses its second zero. Between the two there is the correct value for which the held 
stops on the opposite slope of the barrier as r crosses its second zero. We hnd this correct 
value corresponding to CDL by a simple undershoot/overshoot method. 

After hnding the CDL solution for 0(r) and r(r) we numerically perform the action 
integral in a form simplified using using equations of motion to finally get 

S^cdl] = 47 r 2 J dr (r 3 C(^(r) - 3r(r)) . (5.11) 

Which is the hnal result used in (5.5). 

6 Conclusions 

In this paper we have analysed the evolution of the universe for the Starobinsky potential 
with higher order terms parametrised by two constants Ai and A 2 . In the Sec. 2 the general 
features of the Einstein frame potential have been discussed. Among them the presence of the 
plateau and of the steep slope, the existence of a minimum at (j ) m ; n and of a local maximum 
at (j > ma v (for Ai < — A s ) or a saddle point for Ai = —A s . We have analysed the behaviour 
of the field near extreme points and established possible scenarios of the evolution of the 
system. In the special case of A 2 = 0 the equivalent f(R ) theory which corresponds to the 
model considered in this paper has been found. It has been shown that higher order cor¬ 
rections to the potential weaken the problem of initial conditions of the Starobinsky inflation. 

In the Sec. 3 the inflationary evolution of the system has been outlined together with 
the discussion of the power spectrum of initial inhomogeneities. The issue of eternal inflation 
has been addressed. In fact, the existence of the steep slope has been proven to provide a 
solution to that problem for Ai > 10~ 3 or A 2 > 10 -4 . The attractor solution consistent with 
the slow-roll approximation has been determined. It is separated into two parts: one which 
leads to inflation on the Starobinsky plateau and one which leads to the minimum at cj> m i n . 
Those parts are separated by a local maximum 0 max , around which quantum fluctuations 
of 0 dominate the evolution of the inflaton. This part of an attractor may be a source of 
topological inflation. We have found the maximal amount of e-folds which can be produced 
for general set of A parameters together with a peak of IV max , which corresponds to the saddle 
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point inflation. 


In the Sec. 4 it has been proven that thermal corrections to the Einstein frame poten¬ 
tial shall not modify significantly the evolution of the field. The inflaton trapped in </> m i n 
cannot leave the minimum via thermal correction because non-zero temperature can never 
(even for T beyond the Planck scale) fill the minimum and push 0 towards the GR vacuum. 
Temperatures larger than M p make the minimum even deeper than it is at T = 0. 

In the Sec. 5 the issue of quantum tunnelling between the local minimum at 4 > m ; n and 
the GR vacuum has been analysed. It has been proven that inflationary potentials which 
are able to generate IV max > 60 can generate tunnelling only via the Hawking-Moss effect, 
which gives a lifetime much bigger than the age of the Universe. For Ai 10 6 , which is 
highly inconsistent with inflation, one obtains tunnelling via Coleman - De Luccia instantons 
and oscillating instantons, but the lifetime of the vacuum at cj) m j n is again much bigger than 
the age of the universe. Thus the minimum at is stable and <j) can leave it only via the 
classical evolution. 

Acknowledgments 

This work was supported by the Foundation for Polish Science International PhD Projects 
Programme co-financed by the EU European Regional Development Fund. 

This work was partially supported by National Science Centre, Poland under 

research grants DEC-2012/04/A/ST2/00099 and DEC-2014/13/N/ST2/02712, and grant 

FUGA UMO-2014/12/S/ST2/00243. 


References 

[1] A. A. Starobinsky, Phys. Lett. B 91 (1980) 99. 

[2] P. A. R. Ade et al. [Planck Collaboration], Astron. Astrophys. 571 (2014) A22 
[arXiv:1303.5082 [astro-ph.CO]]. 

[3] A. De Felice and S. Tsujikawa, Living Rev. Rel. 13 (2010) 3 [arXiv:1002.4928 [gr-qc]]. 

[4] A. Codello, J. Joergensen, F. Sannino and O. Svendsen, arXiv: 1404.3558 [hep-ph]. 

[5] I. Ben-Dayan, S. Jing, M. Torabian, A. Westphal and L. Zarate, arXiv:1404.7349 [hep-th]. 

[6] M. Artymowski and Z. Lalak, JCAP09(2014)036 [arXiv:1405.7818 [hep-th]]. 

[7] M. Artymowski, Z. Lalak and M. Lewicki, arXiv:1412.8075 [hep-th]. 

[8] L. Sebastiani, G. Cognola, R. Myrzakulov, S. D. Odintsov and S. Zerbini, Phys. Rev. D 89 
(2014) 023518 [arXiv:1311.0744 [gr-qc]]. 

[9] H. Motohashi, arXiv:1411.2972 [astro-ph.CO]. 

[10] R. Kallosh, A. Linde and D. Roest, Phys. Rev. Lett. 112 (2014) 1, 011303 [arXiv:1310.3950 
[hep-th]]. 

[11] R. Kallosh, A. Linde and D. Roest, JHEP 1311 (2013) 198 [arXiv:1311.0472 [hep-th]]. 

[12] R. Kallosh, A. Linde and D. Roest, JHEP 1408 (2014) 052 [arXiv:1405.3646 [hep-th]]. 

[13] R. Kallosh, A. Linde and D. Roest, JHEP 1409 (2014) 062 [arXiv: 1407.4471 [hep-th]]. 

[14] M. Galante, R. Kallosh, A. Linde and D. Roest, arXiv:1412.3797 [hep-th]. 


18 - 


[15] B. J. Broy, D. Roest and A. Westphal, arXiv: 1408.5904 [hep-th]. 

[16] K. Kamada and J. Yokoyama, Phys. Rev. D 90 (2014) 10, 103520 [arXiv: 1405.6732 [hep-th]]. 

[17] A. H. Guth, J. Phys. A 40 (2007) 6811 [hep-tli/0702178 [HEP-TH]]. 

[18] P. A. R. Ade et al. [Planck Collaboration], arXiv:1303.5076 [astro-ph.CO]. 

[19] P. A. RAde et al. [BICEP2 Collaboration], arXiv: 1403.4302 [astro-ph.CO]. 

[20] M. Bojowald, Living Rev. Rel. 8 (2005) 11 [gr-qc/0601085]. 

[21] M. Artymowski, Y. Ma and X. Zhang, Phys. Rev. D 88 (2013) 10, 104010 [arXiv:1309.3045 
[gr-qc]]. 

[22] X. Zhang, Y. Ma and M. Artymowski, Phys. Rev. D 87 (2013) 8, 084024 [arXiv:1211.4183 
[gr-qc]]. 

[23] A. Ijjas, P. J. Steinhardt and A. Loeb, Phys. Lett. B 723 (2013) 261 [arXiv:1304.2785 
[astro-ph.CO]]. 

[24] L. Kofrnan, A. D. Linde and V. F. Mukhanov, JHEP 0210 (2002) 057 [hep-th/0206088]. 

[25] F. Bezrukov, J. Rubio and M. Shaposhnikov, arXiv:1412.3811 [hep-ph]. 

[26] A. D. Linde, Contemp. Concepts Phys. 5 (1990) 1 [hep-th/0503203]. 

[27] J. I. Kapusta and C. Gale, Cambridge, UK: Univ. Pr. (2006) 428 p 

[28] S. W. Hawking and I. G. Moss, Phys. Lett. B 110 (1982) 35. 

[29] S. R. Coleman and F. De Luccia, Phys. Rev. D 21 (1980) 3305. 

[30] J. C. Hackworth and E. J. Weinberg, Phys. Rev. D 71 (2005) 044014 [hep-th/0410142], 

[31] A. H. Guth and E. J. Weinberg, Nucl. Phys. B 212 (1983) 321. 



